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Abstract. In this paper, we show that astrophysical accretion disks are dynamically unstable 
to non-axisymmetric disturbances. This instability is present in any stably stratified anticy- 
clonically sheared flow as soon as the angular velocity decreases outwards. In the large Froude 
number limit, the maximal growth rate is proportional to the angular rotation velocity, and is 
independent of the stratification. In the low Froude number limit, it decreases like the inverse 
of the Froude number, thereby vanishing for unstratified, centrigugally stable flows. The insta- 
bility is not sensitive to disk boundaries. We discuss the possible significance of our result, and 
its implications on the turbulent state achieved by the disks. We conclude that this linear insta- 
bility is one of the best candidates for the source of turbulence in geometrically thin disks, and 
that magnetic fields can be safely ignored when studying their turbulent state. The relevance 
of the instability for thick disks or nearly neutrally stratified disks remains to be explored. 
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1. Introduction 

The simplest model of an accretion disk is that of an axisymmetric rotating shear flow in 
hydrostatic vertical equilibrium, with a Keplerian velocity law. The hydrostatic state generates 
a vertical stratification. If this stratification is unstable, it leads to turbulence via convective 
instability. When the stratification is stable, it is generally ignored or thought to be unimportant 
in the stability analysis, under the rationale that it can only stabilize the flow. Ignoring the 
stratification makes the accretion disk look like a simple differentially rotating shear flow, 
with an azimuthal keplerian velocity profile V(r) ~ r -1 / 2 . Its linear stability with respect to 
axisymmetric disturbances is governed by the Rayleigh criterion in the inviscid limit: 

for stability. Flows obeying this criterion are called centrifugally stable. The keplerian flow, in 
which angular momentum increases outwards, falls into this category. Yet, there are observa- 
tional evidences that astrophysical (putatively keplerian) disks are turbulent, and thus that a 
source of instability exists in these flows. 

Various mechanism have been found able to destabilize centrifugally stable flows. They may 
or may not apply to astrophysical disks. 

i) Centrifugally stable flows can be destabilized by non-axisymmetric instabilities (Papaloizou 
& Pringle 1984), but these instabilities generally involve reflecting boundary conditions over 
sharp edges, which appear unrealistic in the context of astrophysical disks. 

ii) Centrifugally stable flows can also be destabilized by finite amplitude disturbances involv- 
ing a non-linear mechanism not captured by the Rayleigh criterion (Dubrulle 1993; Richard 

6 Zahn 1999). Laboratory experiments have shown that there is a critical Reynolds number 
above which the flow becomes turbulent (Wendt 1933; Taylor 1936): 

„ rQAr „ /n . 
R = — — > R c . (2) 

Here, Ar is the radial extent of the flow (the gap), Q = V/r is the rotation, and v is the 
viscosity. The critical Reynolds number increases with gap width. For Ar/r — ► 0, it tends to 
that of plane Couette flow: R c ps 2, 000, and for large enough gap (Ar/r > 1/20), the instability 
criterion becomes : 
r 3 Afi 

7 a7 > K (3) 
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this critical gradient Reynolds number R* is of order 6 x 10 5 when the inner cylinder is at rest, 
which is the only case which has been thoroughly explored. 

iii) Centrifugally stable flows can be further destabilized via a physical catalyzer, an ad- 
ditional ingredient making the flow unstable at smaller Reynolds numbers. A first instance 
is a vertical magnetic field (Chandrasekhar 1961). This magnetic field provides a linear axi- 
symmetric instability mechanism for centrifugally stable flows via an interchange instability if 
d{V/r) 2 < 0, i.e. for anticyclonic sheared flows. The application of this mechanism to disks 
was first discussed by Balbus and Hawley (1991); they showed that the stratification of the 
disk does not modify the result, and that the maximal growth rate of instability in that case is 
proportional to the angular rotation velocity in the inviscid limit. The influence of viscosity and 
magnetic diffusivity rj has been recently numerically studied by Riidiger and Zhang (2001); they 
found that instability occurs above a critical Reynolds number, which depends on the magnetic 
Prandtl number P m — v jr\. From their numerical simulations, they fit R c = 51Pm _a65 down 
to Pm = 0.01. In cold disks, Pm = 10 -5 , and so R c is of the order 10 5 (in the inner part of 
the disk) (Rudiger & Zhang 2001). For anyone familiar with MHD flows, the instability of a 
centrifugally stable flow subject to vertical magnetic field is a surprise: magnetic fields generally 
inhibit instabilities, and thus, have a stabilizing influence. In the present case, it appears that 
a stabilizing factor acts upon a stable flow so as to generate instability! 

iv) Very recently, Molemaker et al. (2001) discovered that a similar catalyzing effect could be 
provided by a vertical stable stratification, which induces a linear non-axisymmetric instability 
for all anticyclonically sheared flow. Astrophysical disks are subject to both vertical and radial 
stratification. Our primary goal here is to show that this mechanism may apply to any accretion 
disk and that the inclusion of the radial stratification essentially does not modify the instability. 
Our second goal is to discuss the importance of this mechanism in turbulence generation in 
disks, and compare it to other proposed mechanisms, namely ii) and iii). 



2. A shearing instability in a stratified disk 

2.1. Basic equations 

We consider a rotating compressible stratified disk of finite vertical extent, with velocity 
U = rQ(r)eg and density and pressure po, P- We consider perturbation of the basic state u, p,p. 
In order to eliminate the acoustic waves from the problem, we shall work in the Boussinesq ap- 
proximation. The basic dynamical equations ruling the perturbation are in this approximation: 

V • (pu) = 0, 

Au + u-VU+ — Vp + - — VP = 0, 

Po 5p 

D t h + u • VH = 0. (4) 
Here, D t = dt + U ■ V is the Lagrangian derivative, H = ln(Pp ) and h = — 5p/3po are 
the entropy and the entropy perturbation in the Boussinesq approximation, and assuming for 
simplicity the perfect gas law. Note that the incompressibility condition can be used to eliminate 
the pressure via the Poisson equation: 

— Ap = -V • \D t u + uVU + |^ VP ] . (5) 
Po V 5 Po / 

We work in the local approximation, and decompose the perturbation over Fourier modes with 
wavenumber (k r ,kg,k z ). Because of the mean velocity, the perturbation azimuthal coordinate 
increases with time like: 

9(t) = 9(0) + k e rttt. (6) 
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Because of the shear, the radial wavenumber increases with time according to (Dubrulle & 
Knobloch 1993) 

kr(t) = K(0) - k e St (7) 
where S = rd r Q is the shear rate. To take these two effects into account, we then decompose 
our perturbations along a "travelling" basis exp[ut + irkg(Q — S)t + zk • r], with k • r << 1. 
To leading order in l/(kr) (i.e. considering that S is independent of r), and using cylindrical 
coordinates the Poisson equation becomes: 

Q p 

- ik 2 ^—p = 2VL(k e u - k r v) + s (k r N r N z + k z N 2 z ) , (8) 

where k 2 = k 2 + k 2 + k 2 , u = d z P/d r P(u, v, w) and s = —h/d z H. Here, we have introduced the 
two components of the Brunt- Vaissala frequency: 

N 2 = -—d r Pd r H, 

N 2 = -^-d z Pd z H, (9) 
and used the assumption of rotation on cylinders Q(r) : 

d r Pd z H = d z Pd r H. (10) 
Note that in deriving ([8]) we have taken care to first taking the divergence, and then replace 
the derivative in front of (u, v, w,p) by ik(u, v, w,p). Using @, we obtain the dynamic equation 
under the local approximation as: 

k 2 uu = 2Qv (k 2 - k 2 ) + 2Quk r k e + s (N r N z (k 2 - k 2 r ) + k r k z N 2 ) , 

k 2 uv = -2Qvk r k e + 2Qu(k 2 - k 2 ) - Suk 2 + s (-N r N z k r k e + k e k z N 2 ) , 

k 2 oow = 2Qvk r k z + 2Quk z k e + s (-N r N z k r k z + (k 2 z - k 2 )N 2 ) , 

N r . . 

cos = ——u + w. (11) 

To illustrate the basic mechanism of instability, it is convenient to study first the case where 
N 2 = (no radial stratification). We then return to treat the more general case by building on 
the results of this artificial example. 

2.2. Dispersion relation for N 2 = 



Setting the determinant of flTT| ) to zero, with N r = and N z = N, we find after some algebra 
and rearrangement: 



to 



k 2 + uo 2 (N 2 (k 2 + k 2 ) + n 2 k 2 z + 2nsk 2 ) + 2nSN 2 k 2 = 0. (12) 
Here, k 2 = AQ 2 + 2QS is the epicyclic frequency. In the axi-symmetric limit: k$ — 0, and the 
dispersion relation becomes: 

(k 2 + k 2 z )u 4 + (N 2 k 2 + k 2 z K 2 )uj 2 = 0. (13) 
There are two branches of solution: one with uo = corresponding to neutral modes propagat- 
ing at the mean rotation speed. In the other branch, (k 2 + k 2 )uo 2 = —(k 2 N 2 + k 2 K 2 ). In the 
stably stratified, centrifugally stable case we consider, both iV 2 and k 2 are positive, and u is 
imaginary. This branch corresponds to gravito-inertial waves. As we shall now see, the linear 
non-axisymmetric instability starts from the neutral branch of solution. In the general case (|12D 
is a quadratic equation in oj 2 , which admits two roots, which may be real and distinct, real and 
equal, or complex conjugated. A sufficient condition for instability is therefore that one of the 
root has a positive real part. The root properties are determined by the sign of 2QS = d r Q 2 . If 
this quantity is negative, their product is negative, and the two roots are real of opposite sign, 
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implying instability. In the other case, the two roots are equal or complex conjugated, and their 
real part is negative. A sufficient condition for instability is therefore that 
d r n 2 < 0, (14) 
which is satisfied by Keplerian rotation. The growth rate (to 2 ) 1 / 2 can be found near the neutral 
branch, for vanishing uo 2 : 

2k 2 N 2 ns 

k 2 {N 2 + 2QS) + k 2 K 2 ' y } 

The behavior depends on the Froude number F = Q/N. For strong stratification (F <C 1), the 
growth rate is u r = y/—2Stt, i.e. independent of the stratification. The azimuthal wavenumbers 
in that case scales like k r ~ 1/Ar, while the vertical wavenumber scales like k r /F (Yavneh et 
al. 2001). For large Froude number, the growth rate is 



-2k 2 N 2 ttS 

Since in that case kg ~ k z ~ k r (Yavneh et al. 2001), the growth rate decays like 1/F, in 
agreement with the numerical computations of Molemaker et al. (2001). 

In summary, we have found that for non-axi-symmetric disturbances, the flow is linearly 
unstable, while for axisymmetric disturbances, the flow remains stable. 

2.3. The general case: N 2 ^ 

We consider now the general case. The general dispersion formula is then: 

fcV + u 2 ((k z N r + k r N z f + k 2 z n 2 + k 2 e (N 2 + N 2 + 2QS)) + 2QSN 2 k 2 e = 0. (17) 



In the axisymmetric limit, one obtains a slightly different form of (|i3|): 

(k 2 + k 2 z )u 4 + u 2 ((N r k z + k r N z ) 2 + k 2 z K 2 ) = 0. (18) 
Here we find again a neutral branch of solutions, plus a solution of stable waves in the centrifu- 
gally stable case. 

In the non-axisymmetric limit, the dispersion relation has exactly the same structure as in 
the case with no radial stratification (with the product of the root depending only on SQ), 
and therefore the necessary condition for instability is unchanged. Moreover, in a thin disk disk 
N r ~ (H/r)N z « N z , so that the expression of the growth rates is amost unaffected by the 
radial stratification. 



2.4- Singular limits 

In the present letter, we have used a crude local approximation to derive the stability criterion. 
In fact, these kinds of local approximations have known pitfalls (Dubrulle & Knobloch 1993). 
In the present case, it suggests for example that the growth rate of the perturbation vanishes 
algebraically fast in S and is independent on the stratification. A more careful analysis of the 
normal mode problem however shows that in the small gap limit, and for strong stratification 
(F << 1), the growth rate of the most unstable mode is (Molemaker et al. 2001) 
u w -2S(k e r)e m/s . (19) 
Our analysis and the numerical results of Molemaker et al. (2001) also shows that for large 
values of F, the growth rate behaves linearly with 1/F, thereby vanishing with vanishing 
stratification. This shows that there is no singular limit in this problem. 
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2.5. Non-linear saturation and implications for turbulence 

The present analysis was performed in the inviscid limit. The influence of viscosity and the non- 
linear saturation of the instability have been studied numerically by Molemaker et al. (2001) 
and by Yavneh et al. (2001). They found that a small viscosity essentially does not change the 
results, and only introduces a critical Reynolds number, above which the flow is unstable. For 
example, at F = 0.01, Ro = —2/3, the critical Reynolds number is of order 1200, and from 
Section 2.2 we may extrapolate this result to higher Froude number: R c = 1200V1 + F 2 . In 
thin astrophysical disks, we have N 2 ~ f2 2 (V a d — V), meaning that the Froude number is of 
order 10, if we assume a subadiabatic gradient of 1/10; therefore R c ~ 10 4 . 

The transition to turbulence then proceeds via successive period doubling bifurcations (Mole- 
maker et al. 2001). Without stratification (and without magnetic field), the turbulence would set 
in via finite amplitude disturbances and non-linear mechanisms, a very intermittent scenario, in 
which turbulent domains progressively grow in size with increasing Reynolds number, until they 
finally invade the whole flow, as observed in plane Couette flow (Dauchot & Daviaud 1994). 
Therefore, an interesting open question is what happens in the highly turbulent regime, and 
especially whether this turbulence, triggered by a catalyzer, resembles the turbulence generated 
by non-linear instabilities. 

This question is important, because it is generally believed that turbulent transport properties 
depend on the source of the turbulence. The main argument is that to reach a stationary state, 
turbulence must organize itself so as to suppress the cause of instability. A good example 
is given by centrifugally unstable flows, where the turbulent state is characterized by a flat 
angular momentum distribution (outside boundary layers). But this holds only for moderate 
Reynolds number flow (say up to R = 10 5 ). For larger Reynolds number, there is experimental 
evidence in the rotating shear flow that a transition occurs above which the flow presents all 
the characterisitics of classical shear (boundary layer) turbulence (Lathrop et al. 1992). In 
astrophysical disks, the Reynolds number tends to be much larger than this critical value, and 
one may wonder whether the turbulent transport in disks is not universal. It would therefore 
be very interesting to conduct laboratory and numerically experiments of centrifugally stable, 
stratified flows to study that large Reynolds number turbulent regime. 



3. Conclusion 

We have shown that all accretion disks are subject to a powerful, non-axisymmetric instability. 
This instability lurks in any astrophysical disk, because of the stratification induced by the 
vertical component of gravity. It is a purely hydrodynamical instability and does not require 
the presence of any magnetic field, whatever small. Therefore there is no reason to invoke a 
"dead zone" in insufficiently ionized disks, as was done by Gammie (1996). 

In astrophysical disks, the critical Reynolds number to trigger this instability is of the order of 
10 4 , which is less than the critical Reynolds number for both the magneto-rotational instability 
and the finite-amplitude hydrodynamic instability. But the question of which instability occurs 
first has little relevance: what one should ask is what kind of turbulence develops at high 
Reynolds number. In particular, it is not clear that magnetic field plays an important role in 
that regime, because in most turbulent flows studied so far, the level of magnetic energy turned 
out to be an order of magnitude less than that of kinetic energy. Of course, this does not rule 
out the transport of angular momentum by a large scale magnetic field. 
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